Stripes in high temperature superconductors dilate under 
the influence of an external electric field 
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In this work we study the influence of an electric field on stripes and compare our results with ex- 
perimental findings from scanning tunneling microscope measurements. By introducing a negative- 
bias electric field into a time-dependent Ginzburg-Landau equation with stable stripe solutions, we 
show that hole(electron)-like stripes widen (narrow) compared to the field free case. When a mag- 
netic field is introduced instead, stripe formation is found to be suppressed. 
PACS: 74.20.De, 89.75.Kd, 74.20.-z 



Theoretically stripes in high-temperature supercon- 
ductors (HTSC) have been predicted by Emery and 
Kivelson in Refs. Based on their theoretical 

study on electrons hopping in an antiferromagnetic back- 
ground, Emery and Kivelson found that it is energeti- 
cally favorable for electrons to organize themselves spon- 
taneously into stripes. Later White and Scalapino along 
with Lorenzana and Seibold gave further theoretical sup- 
port to the existence of HTSC stripes in their density- 
matrix renormalization group (DMRG)^ and mean-field 5 
studies. 

Early experimental evidence for the existence of 
stripes came from neutron scattering studies in 
Lai.6- IC Ndo.4Sr £C Cu04 by Tranquada et al. 6 . In their 
later STM measurements of Bi2Sr2CaCu208+<^, Hoff- 
man et al. found density of states (DOS) modulations 
that widen along the (7r, 0) direction and narrow along 
the (7T, 7r) direction, for increasingly negative bias voltage. 
Instead of identifying these DOS modulations as stripes, 
Hoffman et al. concluded that these are quasiparticle in- 
terference patterns, based on comparisons of their STM 
results with ARPES data. However, in a review on fluc- 
tuating stripes in high-temperature superconductors^, 
Kivelson et al. argue that the detection of quasiparti- 
cles does not rule out the existence of stripes. 

In general, STM experimentalists assume implicitly — 
as Hoffman et al. did — that the electric fields generated 
by their STM tip are too weak to influence the surface 
being probed. However, since stripes are ultimately ag- 
gregations of charge carriers 9 , it is much more reason- 
able and intuitive to expect the opposite. In this work, 
we explore how a charged stripe would respond to ex- 
ternal electric fields, and compare our findings with the 
observations in Ref. 0. 

To do so, we can start from a lattice model, for ex- 
ample the Hubbard model or the t-J model with stable 
stripes and then introduce an electric field. However, in 
general it is not possible to solve such models analyti- 
cally. Approximate solutions can be obtained but their 
accuracy cannot be guaranteed, whereas exact numeri- 
cal solutions can only be obtained for small systems at a 
large computational cost. 

An alternative approach relies on the descripion of the 



macroscopic state of a system by an order parameter. 
Such order parameters characterize the thermodynamic 
behaviors of the system of interest, and obey Ginzburg- 
Landau equations (GLE). GLEs without stripe solutions 
have been used to describe superconducting systems^—. 
Of the various equations with stable stripe solutions, the 
complex Ginzburg-Landau equation (CGLE) 

d t V = + (l + ia)V 2 m - (1 + if3)\^\ 2 ^. (1) 

introduced by Newell and Whitehead 13 is one of the most 
well-studied nonlinear equations in the pattern forma- 
tion community. Here a and f3 are phenomenological 
parameters. It has been used to describe phenomena 
from nonlinear waves to second-order phase transitions, 
from superconductivity, superfluidity, and Bose-Einstein 
condensation to liquid crystals and strings in field theory 
(see review by Aranson and Kramer 1 —) . In general GLEs 
have very similar structures and therefore, even though 
the CGLE was not derived from a microscopic model of 
superconductivity, we expect that qualitative conclusions 
derived from the CGLE also apply to other GLEs. 

External electric and magnetic fields are introduced 
into Eq. [T] through minimal coupling which then becomes 
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where qb is the scalar and A is the vector potential. We 
will use the effects of magnetic fields on stripes to demon- 
strate that it is plausible to associate \I/ with the super- 
conducting order parameter. 

We solve Eq. [2] numerically. For this nonlinear par- 
tial differential equation, spectral method a 15 i 16 are the 
most stable and accurate, where they are applicable. In 
these methods, we first Fourier transform the partial dif- 
ferential equation, to get a coupled system of non-linear 
ordinary differential equations of the form u = cu + F(u), 
where c is a constant and F(u) a non-linear term. These 
equations are then integrated using the exponential time 
differencing scheme^ 



u{t n+1 ) = u(t n )e ch + 



r F(u(t n +T))dT, (3) 



where we approximate the non-linear term as 

F = F n + T (F n -F n _ 1 )/h + 0(h 2 ). (4) 

This pseudo-spectral technique 18 works well if we are in- 
terested in field-free solutions or when only an electric 
field is applied. If a vector potential is introduced the dif- 
ferential operator does not reduce to a multiplication in 
Fourier space and thus spectral techniques are not conve- 
nient anymore. To integrate the CGLE in the presence of 
magnetic fields we use a finite-difference approach^ with 
Adam-Bashforth timestepping of third order instead 2 ^. 
In our finite difference approach we take care of the vector 
potential by applying the gauge- invariant ipU method 21 . 
By introducing phase factors of the form 

lF=expj-i^ A x (e,y)de} (5) 
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FIG. 1. When we start from random initial conditions the 
microscopic superconducting pockets act as nucleation seeds. 
These islands of homogenous phase grow up to a critical size 
before they break up into stripes. The eventual width of 
stripes formed depends on the maximum size the islands reach 
before collapsing. 



IX 2/ =exp|— i J A y (x, /i) d/ij- 
we find the simple equality 

(V - iA) 2 = U X d XX (U X ^j) + U y dyy(U y iP) 
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with IX as the conjugate of IX. Because of the discrete 
lattice spacing, it is necessary to define link variables 
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between two neighbouring points. The Laplacian in 
our gauge-invariant finite difference approximation then 
reads 



(9) 



To understand how stripes form in the absence of ex- 
ternal fields we start from random initial conditions for 
^. This describes a superconducting material slightly 
above the critical temperature T c . In this state micro- 
scopic superconducting pockets with finite lifetimes have 
formed everywhere. As the simulation progresses we see 
that these microscopic pockets act as nucleation seeds, 
from which small islands of the homogeneous supercon- 
ducting phase grow (see Fig. [1]). Upon reaching a critical 
size, these islands break up abruptly and evolve into the 
stripe phase. 

The stability of stripes in relation to a phase-separated 
state in which a charge-density wave and a dilute liquid 
state coexist has been analysed in Refs. [22| and l23l for 
spinless fermions in a Hubbard-like model. In this work 
the stripe phase was found to be marginally more sta- 
ble than the phase-separated state at higher hole-doping. 



However, to the best of our knowledge there has been 
no systematic study on the relative stability of the two 
phases in small droplets, and their associated finite size 
effects. In addition to stripes the CGLE also allows the 
formation of hexagonal patterns. A necessary but insuf- 
ficient condition for stable hexagon solutions is for the 
inversion symmetry \I/ =— >• — \I/ to be broke n 24 ! 25 . How- 
ever, we do not expect to see the hexagonal phase in our 
magnetic field simulations if it is slightly favored energet- 
ically, because its symmetry group is incompatible with 
periodic boundary conditions in a rectangular system. If 
it is highly favored energetically, we ought to see multiple 
hexagonal domains within the rectangular system. This 
was not seen in our simulations. 

Stripes formed in our simulations (shown in Fig. [2^) 
are mobile. Since we do not consider anisotropy our 
stripes are different from the highly parallel DOS mod- 
ulations seen in STM experiments. By introducing an 
anisotropic Laplacian or a pinning potential we are able 
to force the formation of parallel stripes which are then 
significantly less mobile than the isotropic ones. Because 
stripes seen in our simulations are not parallel we find the 
typical wavelength by analysing the Fourier spectrum. 
The spectrum of the field free case is shown in Fig. (2}d. 

To investigate the influence of a STM tip on stripes we 
model the STM tip as a charged sphere with a radius of 
20 nm at a distance 1 nm from the superconducting sur- 
face. In the presence of external electric field the initial 
nucleation phase is identical to the field-free case, if we 
start from the same random initial condition. The field 
strength only determines the critical size of the islands. 
For a positively (negatively) charged order parameter we 
found that islands reach smaller (larger) critical sizes un- 
der a negative bias voltage of a few pico-Volt. For a 
positive bias of the same strength the opposite is true. 
The width of stripes then depends on the critical size of 
the collapsing island. In general we found that larger 
(smaller) islands result in wider (more narrow) stripes 
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FIG. 2. (a) Ginzburg-Landau order parameter of the stripe 
solution in the absence of fields and (b) the amplitudes of 
Fourier coefficients. There are no preferred orientations for 
the stripes in the isotropic CGLE, Eq. [2] 

(Fig. [3^ &HK). The change of the Fourier coefficients 
for a hole- and an electron-like order parameter under 
a negative bias voltage can be seen in Fig. [3]o and [4]d 
respectively. 

When we compare our simulations with results de- 
scribed by Hoffman et al. — they see narrowing along 
(7r, 7r) and widening along (7r, 0) — we find that the 
dependence on the direction in /c-space corresponds to 
electron- or hole-like character of our order parameter. 
A variety of experiments analysing the band structure 
of cuprates came to different conclusions. Whether the 
Fermi surface is electon-like and centred at (0, 0)2£ or 
hole-like with its centre at (tt, tt)- 7 is unclear. With our 
simulation we can qualitatively reproduce the change in 
the width of stripes as observed in Ref. [j if we assume 
that the band structure is hole-like along the (tt, 0) di- 
rection and electron-like along the (tt, tt) direction. 




FIG. 3. (a) A hole- like Ginzburg-Landau order parameter 
under the influence of a negative bias voltage forms wider 
stripes, (b) The Fourier amplitudes then become shifted to 
smaller values compared to the field-free case. 

Finally we check whether it is plausible to associate 
^ with the superconducting order parameter in HTSC 
materials. Since strong magnetic fields destroy super- 
conductivity, we ran simulations using the finite differ- 
ence approach over a range of magnetic field strengths. 
When we apply a perpendicular magnetic (of the order 
of pico-Tesla) starting from a random initial condition no 
stripes are formed and the solution of the CGLE becomes 
identically zero. We then turn on the magnetic field af- 




FIG. 4. (a) For an electron-like order parameter we found that 
stripes becomes narrower under the influence of a negative 
bias, (b) The amplitudes of the Fourier coefficients are then 
shifted to larger k- values. 



ter stripes have formed (Fig. [5jt) , and observe that the 
stripes narrow until they break and form a uniform array 
of peaks (Fig. (5}d). With time the peaks shrink until they 
become washed out completely. It is important to em- 
phasize that these are not the vortices usually observed in 
type-II superconductors. We suspect that these features 
might have become the hexagonal phase if we would have 
used the correct boundary conditions. 




FIG. 5. Starting from well developed stripes, we see that a 
magnetic field destroys the stripe phase. Within times that 
are short compared to the typical time scale of stripe dynam- 
ics, we observe the formation of uniformly- spaced vortex- like 
structures across the surface. Both figures show the absolute 
values of the numerical solutions. 

To summarise, in this work we investigated the re- 
sponse of a charged order parameter with stripes to elec- 
tric fields, by exploring numerical solutions of the CGLE. 
For a hole-like order parameter under the influence of a 
negative bias voltage, stripes in our simulation widen. 
The same trend is observed in STM measurment 7 along 
the (tt, 0) direction. With an electron-like order parame- 
ter we find narrower stripes — a behaviour that is actu- 
ally seen along the (tt, tt) direction in Ref. 0. We further 
solve the CGLE in the presence of a magnetic field, to 
show that both the homogeneous and stripe phases are 
destroyed, and thus it is plausible for \I/ to be associ- 
ated with a superconducting order parameter. While we 
cannot claim that the CGLE accurately describes the su- 
perconducting order parameter in HTSC, we believe the 
qualitative response of the true order parameter to an 
external electric field will be similar. 
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